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TECBNICAL NOTE 2056 

TELOCITY DISTRIBOTION ON WING SECTIONS OF ARBITRARY 
SHAPE IN COMPRESSIBLE. POTENTIAL FLOW 


III - CIRCtllATORY FLOWS OBEYING THE SIMPLIFIED 
DENSITY-SPEED RELATION 


By Lipman. Bers 


SUMMARY 


The method of confuting velocity and pressure distributions along 
wing profiles under, the assumption of the simplified density-speed 
relation, outlined in NACA Technical Note 1006, is extended to the case 
of a nonsymmetrical profile and a flow with circulation. The shape of 
the profile j the speed of the undisturbed flow, and a parameter deter- 
mining the angle of attack may be prescribed. The problem is reduced 
to a nonlinear integral equation which can be solved numerically by an 
iteration method. A numerical example is given.’ 


INTRODUCTION 


This paper treats the flow of a compressible fluid past a wing 
section under the assumption of Chaplygin’s simplified density-speed 
relation (references 1, 2, and 3 ). The method is sufficiently well 
known to preclude the necessity of a detailed discussion. It will 
suffice to recall that it consists of replacing the "exact" density- 
speed relation in a potential flow 
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(where p is the density, q the speed, a the speed of sound, 

7 the ratio of specific heats, and the subscript zero refers to the 
stagnation values) by the "approximate" relation 
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P = 



( 2 ) 


which can "be obtained from equation (l) hy setting 7 = —1. Under the 
assumption of relation (2), the continuity equation for the potential 
of a steady two-dimensional flow 
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becomes the classical equation of a minimal surface. 

For flows past airfoils obeying the simplified density— speed 
relation, the "inverse problem" (construction of a flow past a profile 
without predetermining the shape of the profile) was solved by Tsien 
(reference A) for flows without circulation. A formula generating 
circulatory flows was given by the author (reference 5) and, in a more 
elegant and general form, by Gelbart (reference 6). The latter result 
was also obtained independently by Lin (reference 7) • In a recent 
report (reference 8) the "direct problem" (construction of a flow past 
a given profile) was solved for the case of a circulation-free flow and 
a symmetrical profile. The solution for the general case is given in 
this report. It will be seen that the corresponding boundary— value 
problem is equivalent to a mapping problem, similar to the conformal 
mapping problem occurring in the theory of incompressible fluid. This 
mapping problem may be reduced to an integral equation somewhat similar 
to the well-known equation of Theodors en and Garrick but not identical 
with it even in the case of infinitely slow (and therefore incompress- 
ible) flows. The integral equation can be solved numerically by an 
iteration method which seems to converge, though a rigorous convergence 
proof is still lacking. 

The procedure for computing velocity distribution is described in 
the m«-?n part of the paper; the mathematical derivation and justifi- 
cation will be found in the appendix. 

The method described in this report could be extended to the case 
of gases obeying the actual equation of state (or, for that matter, 
any prescribed pressure— density relation) as was done for the case of 
symmetrical flows (reference 9). The computational labor involved in 
treating this case would be so extensive that a detailed description 
does not seem to be called for at present. 
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Tile author wants to draw attention to a paper "by A. Gelbart and. 

3). Resch (reference 10 ) in which a different method of obtaining 
velocity distributions along preassigned profiles is used. While 
Gelbart* s method does not aim at obtaining exact values', it does 
achieve good approximations and involves very little computational 
labor. 

This investigation was conducted at Syracuse University under the 
sponsorship and with the financial assistance of the Rational Advisory 
Committee for Aeronautics. 

The author is greatly indebted to Miss Elizabeth H. Wetherell for 
conpetent assistance. 


SYMBOLS 


A(») 


B(<d) 

c 

ds 

dS 

E(P) 

f(a) 

*n(®> 

f*((J>) 

Si 

8ik 

G 

h(oi) 


auxiliary function defined by equation (A 23 ) 
local speed of sound 
speed of sound at stagnation point 
auxiliary function defined by equation (A 23 ) 
positive constant 

line element in z— plane; line element of profile P 

line element on minimal surface 0 (x,y) whose 
projection is ds 

domain exterior to profile 'P 

function equal to f*(u> — cdq) 

nth approximation to function f(co) 

function defining mapping- of circle into profile P 

function defined in section 3 under .ANALYSIS 

coefficients of metric dS^ 

complex potential of a compressible flow 

function defined by equation ( 10 ) 
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L* 

Moo 

P 

q. 

i«> 

q*[f(«>)] 

Afv 

Q* 

^ ttirt 

a* 

^ 00 

E 

Ee 

s 

S 

U,T 

V 

W* 

z = x + iy 
Z L 

Zfp 

Z(tr) 


line through. Z 5 directed toward E(P) and "bisecting 
angle at 

function defined hy equation (A35) 
stream Mach, number 
profile in z-plane 
speed 

speed of undisturbed flow 
value of q. at a boundary point 
distorted speed 

value of q* at a boundary point 
maximum of *q* 

distorted speed of flow at infinity 
radius of circle in £ —plane 
real part of following term 
arc length measured along P 
total length, of curve P 
components of velocity 
conplex velocity 
distorted complex velocity 
Cartesian coordinates in z-plane 
complex variable 

complex variable on leading edge 
coup lex variable on trailing edge 

coordinate of profile as a function of parameter a 


MCA TU 2056 


5 


a 

P 

7 

r 

5 = i + iT i 
e 

#v 

0 

e(«x) 

X 

A(cd) 

S/l 

P 

\ 

p 0 

a 

ff l 

} 

0(a) 

x(S) 

V 

a> 


“l 


angle of attack 

angle at trailing edge 

exponent in adiabatic relation 

circulation of flow 

auxiliary complex variable 

slope of Telocity vector 

value of 6 on "boundary 

function determining shape of curve P 

square of distorted speed of undisturbed flow 

function defined "by equations ( 9 ) 

Cartesian coordinates in £— plane 
density 

stagnation density 

dimensionless length, parameter along profile P 
parameter value corresponding to point on P 

prescribed value of f (o>j_) 
velocity potential 
value of 0 at "boundary 

auxiliary analytic function defined "by equation (A22) 
stream function 

argument of a point on circle j £ | = E 

parameter determining angle of attack and circulation 

nth approximation to parameter 0 )q 
value of oj for which f(o>) = cr^ 
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ANALYSIS 

1. The Boundary-Value Problem and the Mapping Problem 


Given a profile P (see fig. 1) in the plane of the complex vari- 
able z = x + iy. The profile P is assumed to be a smooth curve 
except perhaps for a sharp comer at the trailing edge Zip. Let l be 
a directed straight line passing through the point zj, pointing toward 
the do main exterior to P and bisecting the angle at zp. (If P has 
a cusp at zp, l shall be tangent to P at this point. ) The posi- 
tion of P can be specified by prescribing the angle a by which l 
mus t be turned in order to mate it coincide with the positive x-axis. 
This angle will be called the angle of attack and it will be assumed 

that | a | < ,jj. 

Let s denote the arc length on P measured in the counter- 
clockwise direction from the point Zp. It will be convenient to use 
the dimensionless parameter 


a = 2ns /S 

where S is the total length of the curve P. The shape, of P is 
determined by the function e(ff),0 < a < 2jt, which denotes the angle 
between l and the tangent to P at a point corresponding to the 
parameter value a, the tangent pointing in the direction of 
increasing a (cf. fig. 1). Note that 

0(0) = « — j| 

0(2n) = 2n + | 

p being the angle at the trailing edge. 

It is required to find a steady compressible potential flow around 
the profile P which is parallel to the x-axis far away from the 
profile and has there a prescribed speed q^. The flow should obey the 
simplified density-speed relation (2) and satisfy the Eutt a-J oukowski 
condition (i.e., the trailing edge should be a stagnation point, or at 
least, if P has a cusp at zp, a point where the streamline divides 
itself into two branches). It is assumed that the values of a 
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and- <1 00 are sucil that there exists one more stagnation point, at some 
point z-q of P corresponding to the parameter value o^. 

It will he shown in the appendix that this problem is equivalent 
to the determination of a mapping of the domain E(P) exterior to the 
profile P onto the domai n | £ | > R in the plane of the auxiliary 
complex variable £, in such a way that the potential and stream func- 
tion of the flow become conjugate harmonic functions in the £ -plane . 

The mapping should take z = » into £ = » and should preserve the 
length and direction of a horizontal line element at infini ty. 

By this mapping the profile P goes over into the circle | £ | = E, 
the points z = z^ and z = z^ being taken into £ = Re~ icu O 

and £ = -fie^O, respectively (see fig. 2). The point— to— point corre- 
spondence between P and the circle is described by an increasing 
function 


a = f*(ai ) , — cdq < o> < 2jr — o> 0 

such that z = Z ( 0 ) corresponds to £ = Ee icu . In particular 

f *(-“*)) = 0 

f*(ic + co 0 ) = (r L • 
f*(2jt - = 2 * 

It will be more convenient to use the function 

f ( 00 ) = f*(o> — (jDq) 




(5) 


( 6 ) 


satisfying the relations 


f(0) = 0 


f(* + 2 = a L ► 
f(2jr) = 2 it 


(7) 


It will be seen that the- functions 
solving a nonlinear integral equation. 


f(o>) can be computed by 
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2. Confutations of the Telocity Distribution, the 
Angle of Attach, and. the Circulation , 

It turns out that the function f(o>) let ermines completely the 
mapping from the £~i plane into the z— plane, as well as the velocity 
distribution of the flow. The "boundary values of quantities character- 
izing the flow will "be denoted "by tildes and considered as functions of 
the par ame ter c. Thus ^(o) denotes the value of the potential at the 
point corresponding to the parameter value cr. 

The speed of the flow at infinity will "be characterized "by means of 
the parameter 


2 


X = 


V 


(8) 


Set 


A (cd) = © [f ((d)] - P to, 0 < (D < 2it 


2n 


A ((d) =A((d + 2ir) 


and 


(9) 


h(oi) = i [A((d + t) - A((d - t)l cot 


dt 


( 10 ) 


The speed of the flow at the profile is given "by the formula 


q[f((D)] = a Q 


2q*[f((p)] 

1 - q*[f((D)] 2 


( 11 ) 


where 


l+r 

5*[f (<»)] - VX 2 


l £ 


s in — 
2 


(D - 2 (D r 


COS 


h((D) 


( 12 ) 
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The angle of attack a is determined by 

Jp2it 

e[f(a>)] due 

0 


( 13 ) 


Finally, the circulation of the flow equals 


2 n n\fx S sin cdq 

r = -a°— — (14) 


02it 

1 . CO 

sin - 

1-1 

it 

■20 

. , % 2+— — 

e-h(o3.)_ x2 * 

. 03 

sm H 

20 

it 2 

cos 

co-2co 0 h((a) 

0 


“ 





. It will he seen that the function f(co) can he compiited. Thus 
it is possible to determine a flow for a given function ©(c), that is, 
for a profile of given shape and for a given value of the parameter 1 , 
that is, for a given speed of the undisturbed flow. The parameter co Q J 

determines the angle of attack and the circulation. In particular, the 
value cOq = 0 leads to a circulation— free flow. 


3- The Integral Equation 
The function f(co) satisfies the relation 


*10) 


f (to) = 2jt- 


10 



n 0 

2 2P 


2£ 

CO* 

sin — 


e -*(o>*)_ x2 + * 

CO* 

sin — 

”0032 

2 


- 

2 

2 


cko’ 


’2 rt 


i-l 


sin 


or 


3 -h(co» ) _ x2 2+ it 


20 


sm 


co' 


20 

* p 0)*-2a) 0 h(co») 

cos^ ^ e 

2 


■(15) 


dco* 


where h(co) is given by equations ( 9 ) and (10). Since the right -han d 
side depends on f, equation ( 15 ) is a nonlinear integral equation. 


For cop - 0 the integral equation reduces to the one derived 
previously (reference 8 ) under the special assumption that the profile 
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is symmetrical with, respect to the x— axis. The special case 1—0 
corresponds to an infinitely slow flow. For 1 = 0 the integral 
equation "becomes 


f(cn) = 2 * 




lo 


CD 1 

sin — 
2 


1 - 


V^® 1 } acD« 


‘2 it 


.1-4 


sin 


< 0 * 


V h( “' > dm> 


10 


It is seen that f(co) does not depend on (Dq. This had to "be 

expected; for in this case the function f (oi) describes the corre- 
spondence "between P and | £ j = E resulting from a conformal mapping 
of E(P) onto | £[ > E. In fact/ an infinitely slow compressible flow 
is equivalent to an incompressible flow. 

The right-hand side of equation ( 15 ) will be denoted 
by f[cDj f(u>* ), X 3 cdq]. The integral equation may be written in the 

symbolic form 

f = F(f, cuq) (16) 


In this equation 1 and cdq are fixed parameters. Instead of deter- 
mining the speed of the flow by the value of 1 (i.e., by the value 

of q ) it may be desirable to prescribe the value of the maximum local 
speed” q max . It is known (and follows from the formulas given by the 
appendix) that the speed attains its maximum at the boundary. According 
to equation (ll), prescribing the maximum local speed is tantamount to 
prescribing the maximum of the function q* defined by equation ( 12 ). 
Let this maximum be denoted by l*juax’ T°ll° ws from equation (12) 


that 


VX 


i* 


max 



- 

• 

JP 


' 



1 a & 

jr 

-cu — 2 cdq 

h(co) 


max 

Stag 


cos 2 

e 



l 



- 



-1 


(17) 


0 < cd < 2n 


x - g i( f ’ “0) 


or 


(18) 
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where g^ denotes the square of the right-hand side of equation (IT). 

If q. *^ is the prescribed value and' X the unknown quantity, then 
equation (17) must he replaced by the system 


f = F(f, X, cdq) 

1 = 8 ^, a> 0 ) I 


(19) 


Similarly, if instead of prescribing the value of cdq it is 
desired to prescribe the angle of attack a, equation (l6) is to be 
replaced by the system 


f = F(f, X, O3 0 ) 

= S2( f > X* °0 


( 20 ) 


where g^(f, X, a,) is obtained by solving equation (13), 


g 2 (f, X, a) = a + 



0[f(aj)] dco'- 


The angle of attack (and hence also a> 0 ) will be determined if the 
value of f(cn) is prescribed at some fixed point cb-.^ 0 < (Bt < 2it. In 

. /i\l O/iw -L 


fact, expanding the terms cos2 


2 cbo 


in equation (15) by means of 


the addition theorem, it is seen that the condition 


f K) = *1 (21) 

is equivalent to a trigonometric equation for the angle cbq. The 0 
coefficients of this equation will depend upon the function f and 
upon X (as well as upon (Bj_). Denote the solution of this trigono- 
metric equation by 

<u 0 = 83 (f* A, cTp.) (22) 


■r 
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If the angle of attack is determined "by means of condition (2l) the 
equations to he solved take the form 


f = E(f, 1, 05 q ) 

«> 0 = g 3 (fj ®i) 


( 23 ) 


Although this way of fixing the angle of attack seems rather artificial, 
it tur ns out to he advantageous for numerical computations. 


Solution of the Integral Equation 


All four forms of the integral equation given in the preceding 
section suggest the application of the method of successive approxi 
mation. To solve equation (l6) for instance, an initial function fo(a>) 
is selected, such that 

V > 0 


f o (0) = 0 


(2*0 


fg(2jt) = 2n 


The operator E(f, 1, od 0 ) transforms this function into another func- 
tion f 1 (<o) which also satisfies conditions (2 4). Now compute success- 
ively the functions 

fl = F(f 0 ) 


f 2 = rffj 

• • • 

Vl “ '( f n) 


If the sequence (f n ) converges and F(f n ) E(lim f n ), then the 
limit f = Hm f n is a solution of equation ( 1 6 ). 
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It is easy to see how this procedure has to he modified for the 
systems (19)., (20), and (23). in the case of system (23), for Instance, 
an initial ‘function f 0 (oo) satisfying conditions (2h) is chosen, and 
the function f2_(o>) and the values cjdq^) are computed so that 

- f l = “ 0 <:L) ) 

“o (l) - «3( f i' «i) 

This .is clearly possible, since the function f Q determines a func- 
tion h(o>) according to equations (9) and (10). With this 
function h(o>) and the given values of X, and a*, the value o^C 1 ) 
is determined so that the condition 

= a l 

he satisfied. Once the value of is found, the computation 

of ^(cjq) can he completed. The function f^co) satisfies condi- 
tions (24), so that the procedure may he continued. The function f n (cu) 
and the numbers are computed by the scheme: 

f n = • F ( f n-l J *■’ “0^ 0 
^ = s 3( f n' ff l) 

The procedure yields the solution of system (23) provided that it con- 
verges, that is, provided that the limits 

f(o)) = lim f n (o>) 
n — > 00 

, • ( n ) 

Qq = lim co n ' ' 

n — » 00 u 

exist. 

The systems (19) and (20) may he treated in a similar way. 
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Attempts to settle the convergence question theoretically have 
failed thus far. On. the basis of numerical computations which were 
carried out, the author believes the iteration method will in general 
diverge when applied to the equation (l 6) and will converge for the case 

of system (23), provided the value of cr is chosen in such a way that 
the point of P corresponding to the value of cr^ is -close to the 

point of maximum curvature of the profile. The convergence might be 
due to the fact that all successive approximations satisfy condi- 
tion ( 21 ). This opinion is substantiated by the fact that previous 
computations (see reference 8) indicate that the iteration method can 
be applied successfully to equation (16) provided the profile is symmet- 
rical and the flow is circulation— free. In fact, for a symmetrical 
profile the condition cdq = 0 is equivalent to a condition of the form 

of equation (21) with o>]_ = Cp = it. 


The computation of the successive approximations is a routine 
matter involving only numerical integration. Note that the integral in 
equation ( 10 ) is a proper Riemann integral, the value of the integral 
at t = 0 being defined as 


lim • 
t — >0 


(a(o) + t) — A(co — t)J cot i 


= 4 A» (to) 


= 4- 0 t [f(o>)]f*(oi) 


n + 3 

2jt 


(26) 


The rapidity of the convergence will depend upon the choice of the 
approximation of order 0 , f 0 (cu) • In general there should be no diffi- 
culty in finding a good initial approximation. The remarks concerning 
this made in reference 8, section 6, apply to the present case with 
obvious modification. 


5. Example 


As an e xamp le, the velocity distribution of a circulatory flow past 
a symmetrical Joukowski profile has been computed. The stream Mach 
number was chosen to be M^ = O.685. Since the stream Mach number is 
connected with the parameter X by the relation 



( 27 ) 
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the value of - X is O.I57. The thickness 
"by setting the usual thickness parameter 
significance of e is seen from fig. 3.) 
determined by choosing the values 

£1^ = Tt 

.. a 1 = 3.1222 

The computation yielded the values 

= 3 ° 27 * 
a = 2° 27* 

The functions ©(c) and ©‘(a) for a Joukowski profile were 
computed by Saltzer (reference 11 ). The numerical values of these 
functions for a symmetrical profile with e = 0.15 are given in 
reference 8. The solution of the Integral equation 
for X = 0 . 157 , cd 0 = 0 ° was chosen as' f 0 (<n). The values of the 
successive approximations f n (co), cjOq^ 11 ) are given in table 1 . The 
functions f n (co) were computed for a) = 0 °, 10 °, . . . , 360°. The 
difference between fg and f^ may be considered as negligible within 
the limits of accuracy of the computations. The resulting velocity 
distribution is given in table 2 and plotted in figure 3. The points 
on the profile are characterized by means of the amplitudes 5 of the 
points of the circle into which they are taken by the standard conformal 
mapping. 

The resulting velocity distribution may be compared with that of 
an incompressible flow = 0) for the same angle of attack, and 

with that obtained by means of the Xarman— Tsien velocity-correction 
formula 


of the profile is determined 
e equal to 0.15. (The 
The angle of attack was 



( 28 ) 


where the subscripts c and i denote compressible and incompressible 
flow, respectively. 
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BpmsrV — During the present investigation, it has come to the 
author’s attention that, in a previous work on this subject (refer- 
ence 8), there are two errors. _ In the formula (6l) on page 2k (cf. 
equation (28) of the present report) the exponent 2 is missing in the 
denominator. In tahle Illb the values in the second, column have heen 
computed incorrectly. The corrected values are given in tahle 3 of the 
present report. 


■ CONCLUDING- REMARKS 


The integral equation used for the numerical computation of 
velocity distributions can also he used as a basis of existence 
theorems for flows obeying the simplified density— speed relation. 
This, however, is a problem of pure mathematical interest. As stated 
in the INTRODUCTION, the extension of the present method to subsonic 
flows obeying the adiabatic pressure-density relation hinges essen- 
tially on problems of computational technique. 


Syracuse University 
Syracuse, N. Y., 


September 1, 19*4-7 
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APPENDIX 

MATHEMATICAL DERIVATION AND JUSTIFICATION OF PEOCEDURE 
1. The Boundary— Value Problem 


This appendix contains the derivation of the results announced in 
sections 1- to 3 under ANALYSIS . According -to the def ini ti on of th 9 
quantities a, S, 0 , © ihe equation of the profile P may he written 
in the form 


— ia 


z = Z(a) = z T + ±- I e ®(t) dtj 0 < o < 2n 


(Al) 


The velocity potential 0 is defined by the relations 

50 


u = a, 


v = a, 


0 5x 

50 

¥ 


(A2) 


where u and v are the components of the velocity vector. As stated 
in the INTRODUCTION , 0 satisfies the minimal surface equation 


1 + 


50\ 2 5 2 0 

/ 3x 2 


50 50 5 2 0 
5x 5y 5x 5y 


1 + 


50\ 2 5 2 0 

W J^2 


0 


(A3) 


It is required to find' a solution of equation (A3) defined in the 
domain E(P) exterior to P and satisfying the conditions 


50 

5x 


a o 


50 


as 


50 

5n 


= 0 


on P 


z — » “I 




(A h) 
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(n denoting the direction normal to P) as well as the Kutta— Joukowski 
condition 



< 00 at 




In general 0 will not he one— valued. The circulation of the flow is 
given hy 


f = 


(u dx + v dy) = a.Q 


d0 


(A5) 


the integration being performed along any simple closed curve 
containing P in its interior. It. is assumed that there exists exactly 
one stagnation point different from z^ = Z^OjJ. 

Set 

w = u — iv = qe (A6) 

The function 9 (slope of the velocity vector) satisfies the conditions 


0-9 0 as z -9 » 


(A7) 


9 


f@ — a — n an the upper hank of P 


© — a — 2jt on the lower hank of P 


(A8) 


Here the upper or lower hank of P denotes the arc of P corresponding 
to the parameter values 0 < 0 < or-^ or cr^ < cr < 2jt, respectively. 

The stream function i|r of the flow is defined hy the relations 


pu _ St 
p 0 a 0 ^ 

> 


P0 a 0 


8t 

Sx 


(A9) 


where p is given hy equation (2). Since t 
streamline, it may he assumed that f = 0 on 


is constant along any 

P. 
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The "distorted complex velocity" is given "by 


w* = q*e 


-id 


(A10) 


vhere q* is the "distorted speed" 


q* = — 
a 0 


(All) 


1 + 1/1 + 


Rote that the parameter X defined hy equation (8) satisfies 




(AL2) 


where q* is the value of q* for the undisturbed flow, and that 
* 00 t * 

equation (All) is equivalent to equation (ll). 

It has been shown (see, for instance, reference 5 ) that the complex 
potential 

G = j + ii (A13) 


is an analytic function of the complex variable w*. 


2. Mapping of the Profile onto a Circle 


In this section the existence and uniqueness of the mapping of the 
domain E(P) described in section 1 under ANALYSIS will be established. 
This transformation 


5 = 6(x,y) 
q = Ti(x,y) 


(All*-) 


must satisfy the following conditions : 


(l) The transformation (Aik) maps E(P) in a one—to—one way onto 
the domain |£[ > R, where E is an appropriately chosen positive 
constant and £ = g + it). The profile P is taken into the 

circle £ = Be^ 0 , 0 51 < 2«. 
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(2) The transformation (Aik) takes the point z = 00 into £ = «> 
and, as z co. 


|L^i 


dx 

dx 


dq 

dy^° 

dy 


-> 0 


c>0 


(A15) 


(3) In the £— plane the complex variables G- and w* are analytic 
f met ions of £ . 


Let ds he a line element in the z— plane , ds = (dx,dy). Let dS 
denote the length of a line element on the minimal surface 0 = 0(x,y) 
whose projection is ds. 


Then 


as 2 = 10 2 + dx 2 + dy 2 

1 p O 

(u dx + v dy) + ds 

a 0 2 


or 


where 


as 


2 = g lx dx 2 + 2g l2 dx dy + g 22 dy £ 


gH = 1 + (l/a Q ) 2 
Sl2 = (l/ a o) 2 cos 
822 = 1 + (i/ a o) 2 


p 

cos^e 


9 sin 9 
cos 2 0 


By virtue of classical theorems it is possible to map the domain E(P) 
in a one— to— one way onto the domain exterior to a circle in the £— plane, 
the mapping being conformal with respect to the Biemannian metric dS. 

It has been shown (reference 5) that any such mapping satisfies condi- 
tion (3). By a linear transformation of the £— plane it can be achieved 
that the mapping satisfy condition (2). The constant E is then 
uniquely determined. 
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According to the definition of the function f*(o)), equation (4) 
tafces the point z = z[f*(a>)J into £ = Re ico . 


3. The Complex Potential in the £- plane 

In the plane G = 0 + ity is an analytic function for j £ j > R. 
The harmonic function i|r is one-valued, and i|r = 0 for 1 1, | = r. The 
harmonic function (f> increases by . r /&q as £ goes once around the 
circle | £ | = B (in the count erclockwise direction). Furthermore, by 
conditions (A4) and (A15), 

G ’(“> = if (Aifi) 


It follows that G has the form (except perhaps for a nonessential 
additive constant) 





r 

a 0 2rri 



(A17) 


so that 


= To, / _ Re ^ 

d£ a 0 l l 


1 + 


Re 


iV 


where 


. -1 r 

a) n = - sm 1 — — — 
0 ^<1. 


(A18) 


The boundary value of the velocity potential is given by 
?[f(“)] = ^loo/aofcos 03 — 0) sin o)q^ 


so that 


#!W] = _4r sm m + gO 


03 — 


03q 


do) 


l 0 


(A19) 


2 


cos 
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Note that the line = 0 in the £ -plane meets the 

circle (£| = R at t, = Re ia3< “ ) and at £ = — Re 1 whereas the 

line = 0 in the z— plane meets the profile P at Zrp and z^. It 
follows that 


f*(-^o) = 0 


f*(* + cu o) = <*L 


(A20) 


as was asserted in section 1 under ANALYSIS. 


4. The Telocity in the £— Plane 

The distorted Telocity w* is analytic in the £— plane. 
Since q ^ 0 in E(P), v*(£ ) j- 0 for 1 > R, so that' 


cs 


log v* = log q* — i0 

is a regular analytic function in |£ j > R. By virtue of equations (A8) 
and (Al8) the Boundary values of the harmonic function 9 are as 
follows : . 


e[f*M] = 


0[f*((D)] — Jt — <X, -CUq < 0) < Jt + (Dq 


0 [f*(co)] — 2jr — a, it + a> 0 < m < 2jt — 


(A21) 


Thus, if f* is assumed known, 9 may he computed in the whole 
domain |£j >R (say, hy Poisson’s integral formula), and so 
can log q* and therefore q. For £ = °° the mean— value theorem for 
harmonic functions yields 


1 

£ = 00 


1 

2 k 


! ' 2«-o)q 

-<Dq 


.©[f*(a>)] dm - a - ^ + cdq 
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Recalling the definition of f (m) (cf. equation ( 6 )), this may he 
written 


i2jt 


^ I ©[f (as)] dm - a - ^ + (Dq 


£ = “ 

But, hy condition (AT), 


e| = 0 

£ = “ 


Formula ( 13 ) for the angle of attack follows from these two relations. 

It is seen from equations (A21) and ( 3 ) that 8 [f*(co)] experiences 
jumps of magnitude (3 and it at m = -o)q and co = n + coq, 

respectively. This implies that q* vanishes at £ = Be 1CU 0 

0 


as 


£ - Be 


-iCDQ 


and at £ = -Re 


ICO, 


'0 


as 


£ + Re 


io)0 


Hence 


xCO = 


w*£ 


l+ ! 


log x(^e 1 ’' JJ ) = A(o>) + iB(co) 


_io)'\ 


(A22) 


({ - + Se i“o) • 

is regular for | >R, continuous for | t, j > R, and everywhere 
different from 0 or 00 . It follows that log x(£) is' regular 
for | £ [ >R and continuous for J £ | > R. Set 


(A23) 


A classical theorem "(see, for instance, reference 12, p. 2k3) yields 
the relation 


A(cd) = — ~ I [b(<c)+ t) — B(co — t)3 cot — + log |x(°°) 
* JO 

Clearly ^ since q*^ = 


(A 2k) 


(A25) 


X(») = /X 



2k 
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It follows from equation (A22 ) that 


A(cd) = log 


q*[f*(cof] 


(A26) 


l+£ 
2 11 


. <0 + ofc 

sm — - — - 


cos 


CD -V CDQ 


Also 


B(cd) = — 0 [f*(cu)] + ^1 + -^cd — & arg ^e icu — e icu °^ — arg (^e ia) + (A27) 


Since 0[f*(cD)j] is given "by equation (A21) '• 

B(cd) = — a*(cd) + Constant 


where 


A*(cd) = ©[V*(cd)J — 
A*(ci> + 2ic) = A*(a>) 


3 + it 

~2 n“ 


CD 


From equations (A26), (A 28), and (A 2k), it follows that 

3 


3 


q*[f*(cD)] = {X 2 * 


CD + 


sin 




it 


CD — CD, 


COS 


'0 


h*(cD) 


(A28) 


(A29) 


(A30) 


where 


1t*(cd) 


I [a*(cd + t) — A*(cd — t)J cot ^ 

Jo 


dt 


(A31) 


Formula (A30) ^is clearly equivalent to equation (12). 


5. The Mapping Function 

The knowledge of the function f(cD) implies the knowledge of, the 
mapping from the domain j£| > R to the z— plane. In fact, using 
equations (2), (A2), and (A9), it is seen that 
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a0 = ^ (is cos 0 + d y sin 0) 


at = -sZ* 


<2= (— dx sin 0 + dy cos 0) 


1 + 


so that 


dz = 


1 iel 


dx + i dy = - e 

* C. 


d0 + i di|r 

— l*(a0 - i dt) 


or 


l/dG 


dz = =• w* dG- 

2 Vw* ) 


(A32) 


where a "bar denotes the conjugate complex quantity. This well-known 
relation is due to Chaplygin (reference 1). In the f; -plane it becomes 



<*«(£) 

w* 


-&»(£)** d£ 


and since z = Zrp corresponds to £ = R e ~ icD 0 = 

£ ~p? “ 

G 5 (£ f ) cL£* 1 

2 / (A33) 

0 ^0 

/ 

Both integrals are path— independent, the integrands being analytic 
functions . 



an 


In particular, if | £ | = R 
arc of the circle |£| = R. 


the integration can be 
Then 


performed along 


a t; = IRe ia) to 



2 6 
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where a* is given by equation (A30) and 9 by equation (A21). Also, 
by equation (A19).> since ^ = 0 on |£| = E, 


.« ie ^ dSff*(m)| 

E dn 

j,. 03 + 05(1 03 03. 

= hie (q^/ao) sin — g cos — — ! 


If these values are introduced into the right— hand side of equation (A33) 
and it is noted that for £ = Ee i<u , z = z[f*(a>)], the following relation 
is obtained: 


Z[f*(03)] = z t + 


2 


03 


■ia 


\fX 


L*(o3 ’ )e 


i©[f*(03*)] 


du3» (A3h) 


i-O3 0 


where 


L*(03) 


sin 


03 + 03Q 




ffc'M - xi*^ 


03 + 03, 


sin 


m 


1 23 

it 


cos^ 


m ~ ^0 e h*(o3) 

(A35) 


Differentiation of equation (A35) with respect to 03 , together with 
equation (Al), yields the relation 

_P 

JL f*(a)) = 2 L*(tl3) (A 36 ) 

2rt i/T 

Integration of this equation from 03 = -03q to 03 = 2it - 03 q, taking 
into account equations (5 ) j leads to. 


E = 


\[x S 


2 *(<W a o) 


2 Jt— 03 q 


L*(03* ) d03* 


(A37) 


u 
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If this expression is introduced, into equation (Al8), it follows that 


r 



2 +£ 

2 S sin cdq 

2jT-<JDq 

L*(o)) da) 
-u) 0 


(A 38 ) 


It is easy to verify that this is equi valent to equation (lh). 


6. The Integral Equation 


If B in equation (A36) is replaced ty its value from equation (A37) 
and the resulting equation is integrated from to a variable upper 

limit a), it is seen that 


10) 


L*(o)» ) do* 


f*(<n) = 2 rt 


-0), 


’0 


2it-o)o 


L*(a)* ) da)* 


-tl)„ 


(A39) 


Since L* is given hy equations (A 29), (A31) , and (A35), this is an 

integral equation for the function f*(o)). Furthermore, 

since f(o) = f*(o) — oig), equation (A39) is identical with equation ( 15 ) • 


It remains to he shown that the solution of the integral equation 
is equivalent to the solution of the original boundary— value problem. 
To verify this, assume that an increasing function f*(a)) satisfying 
equation (A39) (for given values of X and o) Q ) is given. With this 
function define G-(^) by means of equations (A17), (Al8), and (A 37 ). 
From thfe way the integral equation was set up, it follows that there 
exists an analytic function w*(£) regular for |£j > E and satis- 
fying the equation 


W*(co) = fX (AAO) 


w*(se ia) ) = q*[f*(m)]e 
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where q* is given "by equation (A30) and Q "by equation (A 2 l).. With 
these functions, G and w*, form the mapping function, equation (A 33 ). 
This function maps j £) = E in a one-to-one way onto the profile P. 
Hence 


/ 


dz = 0 


if the integration is performed along | £| = E. By Cauchy's theorem 
the same is true for any single closed curve in |£| > K. Hence 
equation ( A33 ) defines a one— valued transformation from the £ —plane 
to the z— plane. Next, 

8| 3| 3£ 3£ 2 w* 2 

^ + = if! 51 + | gr5»\ 

8q 3n 3£ 3£ \2 w* 2 / 


Thus as £ -> «> 


as *0 a\s/X 

dy 


3 | 

8x 


->0 


~>0 


_> isa i/yfx + 
an a 0 2^ 

By equation (8) this is equivalent to 


81 

8x 

8n 


8y 

85 

8n 


— ^ 0 


c > ■ 0 


as £ co. It follows that £ = °° is taken into z 
the Jacobian of the transformation 


■ Finally 
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j = _ d(z,z) 


is equal to 



(A4l) 


Since it was assumed that f** (cu) > 0, it follows from equation (A39) 
that L*(o)) >0., so that by equations (A30) and (A 35 ) |q*| < 1 . By 

the maximum modulus principle, and by equation (A40). this implies 
that |w*| <1 for | £ | > E, so that 


J > 0 for £ > B 


The preceding statements contain the result that equation (A 33 ) 
gives a one-to^ne mapping of |£| > E onto 'E(P), taking |M = B 
into P and £ = « into z = Therefore the functions 0 and w* 
may be considered as functions of x and y, defined in the 
domain E(P). How, by equations (A 33 ) and (AI 3 ) JJ 


where 


I = Re G(£) 
x = Re H(£) 
7 = Ee K(£) 


H - §/ G '(£ ~ ’"H 
e = - 1 J°'Q* + **) « 


(A42) 


(A43) 



(A44) 


Hence 
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By use of a classical theorem "by Weierstrass, it may be concluded from 
this relation that in the z— plane 0 satisfies the minimal surface 
equation (i.e., is the Telocity potential of a compressible flow). 
Next* 


u — iv = 




£ dz d £ dz / 


(A45) 


Using equation (A33) and the expression (A4l) for the Jacobian, 
easily seen that 

_ 2 w* 

dz 1 — j w*!^ G’ 

> 

_ 2 w* 2 w* 

Sz 1 — w* ^ G‘ 

J 


it is 


(A46) 


On the other hand, 

M = i a . 
35 2 

M.ig. 

3C 2 

so that by equations (A^5) to (Ah-7) 


(A47) 


u — iv = Uq 


2w* 

1 — w* 2 


Con^aring the last expression with equations (A10) and (All), it is 
seen that v* is the "distorted Telocity" of the compressible flow 
generated by Wow equation (A4-0) shows that } satisfies the 

required conditions on the profile P and at infinity. 
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TABLE 1.- SUCCESSIVE APPROXIMATIONS TO THE FUNCTION f (co) 


[joukowski profile, e = 0.15, 1^= O.685, X = 0.157, cu^ = jt, CTj_ 


3 - 122 ] 


CD 

f (®) 

(leg) 

Radians) 

fofa) 

fl(®) 

f 2 (m) 

f 3 (tu) 

f 4(“) 

f 5 (m) 

w 

f 7 ( a >) 

f8(“) 


0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

10 

.175 

.032 

.031 

.034 

.032 

.032 

.033 

.032 

.032 

.032 

.032 

20 

-349 

.128 

.126 

.134 

.128 

.127 

.129 

.127 

.127 

.127 

.127 

30 

.524 

.279 

.277 

•293 

• 279 

.279 

.281 

.279 

• 279 

.279 

.279 

4o 

.698 

.476 

.474 

.501 

.477 

.477 

• 479 

.477 

.477 

.477 

..477 

50 

.873 

.708 

.706 

.745 

.708 

.710 

-711 

• 710 

.710 

.710 

• 710 

60 

1.047 

.962 

.960 

1.011 

.960 

.965 

.965 

.965 

.96 5 

.965 

.965 

70 

1.222 

1.227 

1.225 

1.287 

1.221 

1.232 

1.230 

1.231 

1.231 

1.231 

1.231 

80 

1.396 

1.493 

1.491 

I.562 

1.482 

1.499 

1-495 

1.498 

1.497 

1.497 

1.497 

90 

1.571 

1.751 

1.749 

1.826 

1.734 

1.758 

1.751 

li756 

1.754 

1.755 

1-755 

100 

1.7^5 

1.994 

1.992 

•2.072 

1.971 

2.010 

1.991 

1-999 

1.996 

1-997 

1.997 

110 

1.920 

2.218 

2.216 

2.295 

2.189 

2.242 

2.211 

2.223 

2.218 

2.220 

2.219 

120 

2.094 

2.420 

2.418 

2.491 

2.385 

2.442 

2.408 

2.423 

2.416 

2.419 

2.418 

130 

2.269 

2.597 

2.595 

2.658 

2.558 

2.618 

2.580 

2.599 

2.590 

2.594 

2.592 

140 

2.444 

2.750 

2.743 

2.797 

2.707 

2.768 

2*727 

2.749 

2.738 

2.744 

2.741 

150 

2.6l8 

2.878 

2.867 

2.910 

2.835' 

2.890 

2.851 

2.874 

2.862 

2.869 

2.865 

160 

2.793 

2.983 

2.968 

2.997 

2.945 

2.987 

2.955 

2.976 

2.964 

2.971 

2.967 

170 

2.967 

3.O69 

3.051 

3-064 

3.040 

3.061 

3-044 

3.056, 

3.048 

3.053 

3.050 

180 

3-142 

3.142 

3.122 

3-122 

3.122 

3-122 

3.1-22 

3-122 

3-122 

3-122 

3-122 

190 

3.316 

3-215 

3.195 

3.191 

3.198 

3.191 

3.198 

3.192 

3.196 

3-193 

3.195 

200 

3.491 

3.300 

3.282 

3.286 

3.278 

3-280 

3-282 

3.278 

3.282 

3-279 

3.281 

210 

3.665 

3.405 

3.388 

3.409 

3.376 

3.392 

3.386 

3.387 

3.388 

3.387 

3.387 

220 

3.840 

3.534 

3.515 

3.56O 

3.496 

3.527 

3-514 

3.519 

3.517 

3.518 

3.518 

230 

4.014 

3.686 

3.668 

. 3-737 

3-642 

3.685 

3.668 

3.675 

3.672 

3.674 

3.673 

240 

4.189 

3.863 

3.849 

3-939 

3.819 

3-866 

3-848 

3.856 

3.852 

3.854 

3.853 

250 

4.363 

4.065 

4.055 

4.166 

4.023 

4.-070 

4.052 

4.060 

4.056 

4.058 

4.057 

260 

4.538 

4.289 

4.282 

4.415 

4.251 

4.295 

4.278 

4.286 

4.282 

4.285 

4.283 

270 

4.712 

4.532 

4.528 

4.682 

4.500 

4.538 

4.524 

4.531 

4.527 

4.530 

4.528 

280 

4.887 

4.790 

4.787 

4.962 

4.765 

4.796 

4.784 

. 4.790 

4.786 

4.789 

4.787 

290 

5.062 

5.O56 

5^055 

5-249 

5-039 

5.061 

5.052 

5.056 

5.053 

5.055 

5.054 

300 

5.236 

5.321 

5.321 

5-533 

5-311' 

5.325 

5.318 

5.321 

5.319 

5.321 

5.320 

310 

5.411 

5.575 

5.576 

5.669 

5.571 

5-578 

5-574 

5.576 

5.574 

5.576 

5-575 

320 

5.585 

5.807 

5.808 

5.779 

5.805 

5.809 

5.806 

5.807 

5.807 

5.808 

5.807 

330 

5.760 

6.005 

6.006 

5.989 

6.003 

6.006 

6.004 

6.005 

6. 004 

6.005 

6.005 

340 

6.934 

6.156 

6.158 

6.149 

6.156 

6.157 

6.156 

6.156 

6.156 

6.156 

6.156 

350 

6.109 

6.251 

6.252 

6.249 

6.251 

6.251 

6.251 

6.251 

6.251 

6.251 

6.251 

360 

6.283 

6.283 

6.283 

6.283 

6.283 

6.283 

6.283 

6.283 

6.283 

6.283 

6.283 

“0 


0° 30* 

5° 52 * 

1° 37’ 

4° 47* 
— 

2° 42* 

3° 52' 

3° 17* 

3° 37’ 

3° 27* 
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TABLE 2 .- "VELOCITY DISTRIBUTION 


{joukowski profile, € = 0.15, M*, = O. 685 , 1 = 0.157, = it, a ± = 3 . 122] 


s 

By present 
method. 

For incompressible 
fluid 

By Earman— Tsien 
method. 

0 

0.884 

O.869 

0.831 

10 

.849 

.876 

.840 

20 

.856 

.893 

.860 

30 

.888 

.919 

• 893 

40 

.932 

.952 

• 935 

50 

.984 

• 993 

• 990 

60 

' 1.047 

1.040 

1. 056 

70 

1.121 

1.092 

1.132 

80 

1.204 

1.148 

1.220 

90 

1.295 

I.205 

1.316 

100 

1.392 

1.263 

1.421 

no 

1.492 

1.319 

1.529 

120 

1.589 

1.372 

1.641 

130 

1.672 

1.415 

1.739 

140 

1.715 

1.445 

1.812 

150 

1.675 

1.445 

1.812 

160 

1.464 

1.373 

1.644 

170 

• 994 

1.098 

1.141 

180 . 

•457 

.070 

.059 

190 

.299 

• 377 

.325 

200 

.760 

.837 

• 793 

210 

1.034 

1.047 

1.066 

220 

1.180 

1.141 

1.209 

230 

1.244 

1.177 

1.268 

240 

I.255. 

1.182 

1.277 

250 

1.235 

1.167 

I.251 

260 

1.194 

1.141 

1.208 

270 

1.143 

1.106 

1.155 

280 

1.088 

1.068 

1.097 

290 

1.033 

1.028 

1.039 

300 

.981 

•990 

.986 

310 

• 935 

• 954 

.938 

320 

.895 

.923 

.898 

330 

.863 

.898 

.866 

340 

.839 

.880 

.844 

350 

• 839 

.870 

.832 

360 

.884 

.869 

.831 
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Figure 3.- Velocity distribution along a Joukowski profile, e = 0.15, a = 2° 27', M*, = 0.685. 






